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The shift of the p-wave D s meson mass due to coupling to the DK channel is calculated without 
fitting parameters using the chiral Lagrangian. As a result the original Qq mass 2.490 MeV generi- 
cally calculated in the relativistic quark models is shifted down to the experimental value 2317 MeV. 
With the same Lagrangian the shift of the radial excited 1~ level is much smaller, while the total 
' width r > 100 MeV and the width ratio is in contradiction with the D*(2632) state observed by 

Q ' SELEX group. 

o ' 



Or 

o 



o 
o 



PACS: 11.30.Rd,11.80.Gw,12.38Lg,14.40.Lb 



I. INTRODUCTION 



The heavy-light D and D s mesons have extensively been investigated experimentally in the last 20 years pj. 
Recently a lot of attention has been paid to the + state of the D s meson. It was found by the BABAR group at 
the mass M = 2318± 1.3 MeV and confirmed by CLEO and Belle 0. This value of mass is ~ 40 MeV below the KD 
t— I , threshold and the width is very small. On the theoretical side this state was the object of intensive study (see Ref. Q 
• for reviews and references). There exists some disparity between the theoretical predictions of the D*(0 + ) mass in 
the Relativistic Quark Model (RQM) calculations Q-Pij an d experimental results. Indeed the RQM predictions vary 
from 2380 MeV 0, to 2487 MeV 01 j being however substantially larger than measured experimentally. 

Most recently a new D s state was observed in the SELEX experiment ^l| at the mass value M — 2632 MeV. It is 
argued in Ref. ^2 that this state may be associated with the radial excited 1 _ * level of D s meson, which is shifted 
down due to coupling to the r\D s and KD channels. As one will see both levels D*(2317) and D*(2632) are connected 
to thresholds via chiral decays and we shall treat them below using the chiral Lagrangian containing only f n , fx, fr t 



CL' 

"i ■ as parameters. 

It is the purpose of the present paper to study the effect of nearby thresholds on the position of resonances first in the 
most general setting and then to calculate numerically using the chiral quark Lagrangian without fitting parameters. 
The paper is organized as follows. In section 2 we present a general discussion of channel coupling and level shift using 
£> ■ the (relativistic) Hamiltonian formalism, where we also give a classification of possible 5-matrix poles. In section 3 
the chiral quark Lagrangian is written down and used to describe the decay transition D* — > DK, D s rj. The final 
equation for the resonance position with account of this decay is explicitly written. In section 4 the numerical solution 
of this equation is described and the final results are presented. The paper closes with a discussion and comparison 
with other results. 



II. RESONANCE STATES IN THE COUPLED-CHANNEL SYSTEM 



The relativistic quark model has been remarkably successful in predicting the D and D s meson spectrum, apart 
from some exceptions of a few resonances, which are experimentally found at substantially lower masses. Similar 
results have been found recently in the relativistic Hamiltonian approach^! derived on the basis of Field Correlator 
Method (FCM) and applied to the D,D S mesons in Ref. 15]. The results for the masses in the FCM analysis 
are shown in Table 1 and Table 2. (The entries given in the tables are recalculated for parameters given in the Table 
captions). 
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Table 1 Masses of L = 0,1 states of D mesons. Input parameters used in the 
FCM calculations: a s = 0.46, a = 0.17 GeV 2 , m c = 1.44. GeV, m n = 7 
MeV 



State J p 


0" 


1" 


0+ 


i+(0 




2+ 


1~* 


Mass (MeV) 
from [15] 


1859 


2047 


2370 


2425 


2455 


2456 


2729 


Mass (MeV) 
experiment 


1869 


2010 


2300±60 


2400 


2422 


2459 


(2640?) 


T (MeV) 
experiment 




< 0.13 


280 


- 250 


20 


23H-45 


< 15 



As is seen from Tables 1 and 2 the overall agreement is reasonably good except for a few states. In particular, the 
D* (0 + ) state is one example of such a discrepancy in the prediction of theoretical models and which can be associated 
with the KD threshold at 2366 MeV. 



Table 2 Masses of L = 0, 1 states of D s mesons. Input parameter used in the 
FCM calculations: a s = 0.46, a = 0.17 GeV 2 , m c = 1.44. GeV, m s = 
0.175 GeV 



State J p 


0" 


1" 


0+ 


i+(0 


l+(h) 


2+ 


1"* 


Mass (MeV) 
from [15] 


1929 


2087 


2404 


2462 


2488 


2494 


2774 


Mass (MeV) 
experiment 


1968 


2112 


2317 


2462 


2536 


2572 


2632(?) 


T (MeV) 
experiment 




< 1.9 


< 10 


< 6.6 
8.6 -j- 0.4 


< 2.3 


- 15 


< 17 



For comparison in Table 3 a summary of the results is given of other theoretical quark model predictions for this 
state. A look at the table tells us that all the theoretical predictions are about ~ 90 4- 190 MeV higher. So one needs 
a shift of about this value to get agreement with the experimental value. A similar descrepancy can be seen from 
Table 2 for the D*(2632) resonance. We explore in this paper whether this disparity can be explained due to the 
presence of coupled channels with nearby thresholds. 



Table 3 Theoretical predictions of the mass D* (0 + ) in various quark models 



Ref. 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


Mass (MeV) 


2480 


2388 


2380 


2508 


2455 


2487 



Resonances in the coupled-channel system were considered in numerous papers both in nonrclativistic nuclear 
physics and in the relativistic Hamiltonian dynamics, see [l6| for a review and references. Assuming that a local or 
nonlocal relativistic Hamiltonian can be written for each channel Hi, i = 1,2, ... and for the Channel Coupling (CC), 
Vij, i,j — 1,2, ... the time- independent system of equations can be written as 



[(ff, - E)S W + V u <]Gvr = 1. 



(1) 



3 



For two channels it is 

(H 1 -E)G n + V 12 G 21 = l, (H 1 -E)G 12 + V 12 G 22 =0 1 
(H 2 - E)G 22 + V 21 G 12 = 1, V 21 G n + (H 2 - E)G 21 = 0. 

The system (J3J) can be reduced to the effective one-channel problem, corresponding to the Feshbach equation [l7T | 

{Hx - E)Gu - V12 „ 1 „ V 31 G n = 1. (3) 
ti 2 — tj 

At this point one can classify all possible poles E of the Green functions G^. These poles may originate from the 

(n) 

bound states or resonances in a given channel i, located at El and shifted due to CC to a new position, which we 

will denote by E\ n '* . Another possibility is that resonance poles appear solely due to the strong CC interaction - the 
so-called CC poles 0,0]. These extra poles usually originate from distant dynamical poles in the complex plane, 
which move close to threshold when the CC coupling increases. The quantitative characteristics of the CC interaction 
is given by the last term on the l.h.s. of Eq. ®, which can be called the Feshbach potential, 

Vi2i{E) = -V 12 G 2 V 21 = -Va „ 1 = V 21 . (4) 

ti% — hj 

Note that V\ 2 \(E) can support bound states or resonances even in the case when diagonal interaction Vi, i = 1,2 
vanishes but Vi 2 — V 2 \ is large enough. 

Of special importance for us is the case when in one channel, e.g. i=l, the spectrum is discrete (see Ref. 0] for 
a more extensive discussion), and one is interested in the shift of the discrete level due to the coupling to channel 2, 
where states can be unconfined. 

A somewhat similar approach was undertaken in recent papers jlflj . where in our notations the scattering channel 
2 and the corresponding Feshbach potential V 2 \ 2 was modelled to calculate the scattering cross section in channel 2. 
We shall compare the results of Ref. [l9| with ours in the concluding section. Eq. © connects in general all states 
in channel 1 and channel 2. If one separates one state and neglects all other states in channel 1, then one gets the 
following equation for the position of the pole(s) in the Green function 

E = E{ n) -(n\V 121I -^V 21 \n), (5) 

where E± is the selected unperturbed level in the channel 1. Insertion of the complete set of states \m)(m\ with 
eigenvalues E^"" in the channel 2, yields 

E = E[ n) £<„|tr 12 | m >-l (m\V 21 \n). (6) 

m E 2 — E 

In what follows we shall be using Eq. |J?jJ to calculate the shift AE n — E±* — E± of the cs levels due to the open 
channel 2: KD or r/D s scattering states, neglecting interaction in these states. The most important point is how to 
find the operators Vi 2 . In the next section we shall use the chiral Lagrangian which will provide V\ 2 explicitly without 
free parameters. 

III. COUPLED CHANNELS AND CHIRAL DECAYS 

One starts with the Lagrangian for the flavor SU(3) triplet of quarks in the field of the heavy (c or b) quark [20,21]. 
In the Euclidean notations 

L = i J d 4 .xV + (d + m + M)ip, (7) 

where the mass operator is 

M = mil = M exp (i l5 ^^\ (8) 
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and t a = ^A a , A a is the Gell-Mann matrix, a = 1, ...8, f n — 0.093 GeV and the matrix Nambu-Goldstone SU(3) wave 
function is 



a \a = \/2 



V 



V6 n s/2> 



+ 



K+ \ 



(9) 



M is the (nonlocal) effective mass operator, which in the local limit has the form (see [22,23] for discussion and 
derivation) 



M = crlrl 



(10) 



where |r| is the distance from the light quark (u, d, s,) to the heavy quark (c or b). Thus the Lagrangian (7) contains 
effects of both confinement and chiral symmetry breaking. 

From Eqs. Q and (jHJ), expanding the exponent in Eq. (JSJ), one can derive the meson emission part of Lagrangian, 



AL = — f ^ + (x)a|x| 7 5^^^ 4 x. 



(11) 



This Lagrangian can be expressed as in Ref. [20] in terms of the standard Weinberg Lagrangian [24] . It was used in 
Ref. [25] to calculate the decay widths of heavy-light mesons with good accuracy. 

It is clear that the Lagrangian Q generates (due to the various Fock components in ip + or ip) in general a many- 
channel system of equations for the Green functions. It contains the main channel (e.g. the D s channel) and in 
addition the channel(s) for its virtual decay products like the (D + K) channel or (D s + rj) channel. 

In what follows we shall be working with Eq. © to apply it first of all to the D*(2317) state. In this case 
refers to the + level of the D* system, and E^ refers to the (continuous) energy of the system D + K in the 
orbital S state. One can neglect the DK interaction in the first approximation and write for the wave functions Dirac 
equations 



1 / r (1) o (1) 



(12) 



-Qi -o. 



"jlM 



OAfi 



\m) = 9(D) 



ipr 



(13) 
(14) 



9(D) - ( G(2) ^ 



n( 2 ) _ o q( 2 ) _ o , 



Therefore the matrix elements in Eq. 10 are 



\V 12 \m) = - j 9 + (D* s )a\r\ l5 ^9(D) 



=4 pr 



= d 3 7 



(15) 
(16) 

(17) 



V2a f d 3 r (1 



if* 



Cr7( 1 )+p(2)o+ n, _ p( 1 )+^(2)q+ Qi 



(18) 



\/2ct<5 MiJ 



m 2 / sin pr 



Py / 2e p V 3 
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Now G«,F« and G< 2 ',i?( 2 ) are solutions of the Dirac equation 

^# - fF^ + ( e « _ 7 c (r) - m,)G« - MG« = 0, 
^ + fG<^ - ( E W _ Vc (r) + m,)^" - MF« = 0. 

Here «i = 1, Ki = —1, M = err, V c (r) = — ^r, and 



(19) 



3r 

mi=m s = 0.15 0.25 GeV, m 2 = 0. 



The connection between £n and Ei is 



B{ n) - 4 X) + - C ; 4 m) = e« + m K + m c + = (20) 

If one neglects higher states of the D meson, the sum in Eq.© can be rewritten as 

E = E [n)_y : ^P MVMm)\\ (21) 

1 ^ (2tt) 3 E i { n \p)-E 

From Eqs. (|17|i and l|21l) it is clear that the free Green function of the KD system has the form (we take into account 
the fact that only S- waves of KD are involved) 



G (k, x, aO = (7 ,£? a exp ^r x ' )] * ) 



S— waves 



(22) 



1 rh K sin(fcx< ) cxp(ikx > ) 
47r uj(k) kxx' 



where 



fc 2 = 2m K (E - m D - m K ), u{k) = Jk 2 +m 2 K . (23) 
Finally Eq. (|21|l can be rewritten as 

B = E[ n) - J V(x)V + (x')d 3 xd 3 x'G (k, x, x), (24) 

where V(x) is 

V{x)=Ux)^j^^Mx)= 

(25) 



^( G (i) +i r(2) _ F (i)+ G (2) )n + oMifij 



and we have used relation fli 1M = <xnf2 i 0M . As a first approximation one can use the fact that functions 
are concentrated around the middle point x = b and write 

V(x)^CS(x-b); C= [v(x)dx (26) 
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/ a;y(a;)(ia; 
/y(a;)da; ' 



As a result one obtains 



E = E (n)^ c 2^m K sin(kb)e^^ 
uj(k)k 



where k = ^j2fa^{E — Mp — mx) is the relative momentum of the K meson. 

In the vicinity of the DK threshold one can replace cu(k) = nik (this is implied by the form of Eq. J22J0- Eq. I|27|) 
is a transcendental equation for the position of the pole E. Since C does not depend on k Eq. I|27|l has a simple 
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0.4 



0.3 



0.2 



0.1 



0.35 



0.30 , 




FIG. 1: The trajectory of the pole in the complex fc-plane parameterized by C. For b we have used b = 0.4 GeV . 



square-root threshold at E — mo + mx- The starting position of E = E po i e is at E{ . When one takes into account 
the second term in Eq. (|27|) with gradually increasing C 2 the pole moves to the final value in the upper physical 
k-sheet. One expects that the trajectory will go down in mass, possibly nearby the final value of m{D* s ) = 2317 
MeV. In Fig. 1 we display the trajectory of the pole solution of Eq. (27) in fc-plane parameterized by C. One can 
see that for strong enough channel coupling the resonance pole moves down under the DK threshold, which will be 
substantiated by the exact calculation of V\ i . given by Eq. (18). 

We now turn to the case of the D* (2632) |ll| assuming after Ref. ^2 that it can be associated with the radially 
excited D* s (l~) state. From Table 2 one can see that the expected shift should be around 100-150 MeV downwards, 
and from the channels DK, D s rj the decay is in the p-wave. As before we shall use Eqs. (12) and (17) where now 
instead of Z?*(0 + ) one should write the D*(l~) state, i.e. 



o (1) -n. 



(28) 



The wave function in channel 2 is either the same as in Eqs. 114(1 - l|16JI . or in the case of the 77 channel, one should 
replace D by D Sl and the K meson with momentum p by the rj meson with momentum p'. With these assignments 
of the various states Eq. I|17H retains its form, but Eq. l|18f) becomes 



V2a 



V 3 



(29) 



Upy/2e p V 3 Jo 



, . sinpr 

ar I cos j 

pr 



pr){GW+F<»-FW+GM), 



where p is the momentum of the K or rj meson. 

Here the radial quantum number in Eq. 129|) is the first radial excited state n = 1. To compute the transition 
potential matrix elements for the D* s «-> DK coupling one can use the coordinate representation as in Eqs. i|18[l - H25|) 
including Eq. 1)28(1 . or directly calculate the d 3 p integral in Eq. (|21|l . In the latter way one needs to compute Eq. I|21|l . 

uj{p) = e p = ^p 2 + m 2 K 



E = E 



(n) 



and ui2(p) is 



(2^)3 ,/ 2tu(p)(E 2 (p) - E) 



(30) 



(31) 
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Eq. (30) has to be extended to also include the D .^-channel contribution. Observation of Eqs. © and <|11|) shows 
that in Eq. @ only the lowest diagonal term enters in the transition potential matrix element. Hence we get for 

(n|Vi2|m) an additional factor As a result we have 

The modification of V12, Eq. (29), for the case of the 77-channel is straightforward. Clearly, the assignment of the D* 
state in Eq. Ij28(l remains the same and we have to replace in Eq. (121 H ttlk — > and tod — > mu s . Moreover, the 
state *(JD) in Eq. (17) has to be replaced by ^l(D s ), where F (2) , G (2) now refer to the D s state. 



IV. MASS SHIFT IN THE CHIRAL LAGRANGIAN FORMALISM 



We may determine the energy shift AE using the Dirac wave function as found by solving Eq. (|19|l for the case of 
an effective quark mass operator in the field correlator method |23j |. In the leading order it has a nonlocal form and 
can be parameterized as 



M(x, y) 



1 



2T 3 0F 



exp 



(*-y) s 

4T 2 

a 



(32) 



where a is the string tension and T g the gluon correlation length, characterizing the scale of nonlocality. Note that for 
T g — ► one obtains from Eq. (28) the local limit (10). The physical value of T g found on the lattice and analytically 
is small, T g = 0.25 fm [12]. 

For the case of a single coupled channel D* <-> DK we may write Eq. (30) as 



AF = -J- ( dZ MV(p)\ 2 m , 
(2tt)3 J 2oj(p)(E(p)-E -AEY 1 ' 

2 

where AE — E — Eq with Eq — mp* — mn — rriK is the mass shift and E(p) the kinetic operator E(p) — ■ 

In the considered QCD string model we have taken T g — 0.25 fm, in accordance with lattice gauge simulations 
[T3 |. For a given string tension a the wave functions of the D and D s system can be found as solutions of the 
Schwinger-Dyson-Dirac eauation|23j for the light-heavy quark system. It is given by Eq. (19) with the quark mass 
operator (|32|l . The latter is found from the selfconsistent solution of nonlinear equations, (Eqs. (15), (16) in [22]). It 
exhibits the property of both confinement and chiral symmetry breaking. The states are in general characterized by 
the quantum numbers j, I, k. In particular, the D and the orbitally excited state D* correspond to the solution of the 
ground state in the j,l,K= 1/2, 0,-1 and 1/2, 1, +1 channel respectively. 

In this study a value of a = 0.18 GcV 2 and a s — 0.35 is adopted. For convenience, a zero mass is used for the 
u, d quark, while for the s-quark we have taken m s = 200 MeV. Having constructed these wave functions with these 
parameters we determine the matrix elements . From this we may then solve the resulting Eq. (|33J) iteratively. 
To determine the actual position of the pole we in general have to analytic continue the integral into the second 
sheet in the case that E is above the i^D-threshold. Although this can be done, we will assume in this study that 
the imaginary part of the pole position does not affect the solution substantially , which is certainly true when the 
pole dives under the £) if -threshold, and was checked in other situations. Confining oursclf to real values of E and 
taking the principal value of the integral in Eq. (33) when E is above threshold the energy shift is determined as a 
solution of the resulting equation. In case of a solution above threshold, the width of the resonance can be obtained 
by calculating the discontinuity of the integral at this energy. 

In the calculations we have used for the threshold mass of the KD system tuk + mo = 2.366 GeV. A typical value 
of 2.49 GeV is adopted for the unperturbed D* meson mass. As is seen in Table 3 this is in accordance with the 
predictions of FCM and of many quark constituent quark models. 

In Fig. 2 is shown the prediction of the shifted mass of the D* meson due to the channel coupling to the KD- 
system as a function of the unperturbed mass D s . The flavoured symmetric and broken value of f v — 93 MeV and 
Jk = 121 MeV have been used [1]. The shifted mass is found using 

M Dl = Mqjm + AE, (34) 

where AE satisfies Eq. I|33|) . Due to the coupling to the KD system the D* meson can either become unstable or 
stable, depending on the sign and magnitude of the mass shift. We find in our model, that it can be substantial and 
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of the order of hundred MeV. As is seen from the figure we find the experimentally observed value of 2.317 GeV for 
an unperturbed mass of approximately 2.49 GeV. 

In the string model studied we find that the position of the pole moves well down below the KD threshold, yielding a 
stable state. The channel coupling is found to yield attraction, so that the position of the mass pole shifts downwards. 

Similarly as for the Z?*(2317) meson we may estimate the mass shift of the D*(2632) due to channel coupling. 
The calculations proceeds in the same way. The bare Z?*(2632) is assumed to correspond to the first radial excited 
state with j,l,K = 1/2,0,-1 and to have a mass of 2.76 GeV as found in the FCM Clearly a mass shift of 

about —140 MeV is needed to get agreement with the observed mass. As discussed, coupling can occur in this case 
to KD and rjD s channels. For the D meson we assume the experimental observed mass of 1.869 GeV. There are 
two candidates for the D s state, being the 0~ and 1~ states with masses of 1.968 and 2.120 GeV respectively. Using 
the various wave functions obtained from the QCD string model, the interaction matrix elements are calculated from 
Eq. I|29|) . The two D s states are degenerate in the considered model and as a result have the same wave function, but 
differ in the kinematics of the momenta in view of the adopted mass difference, which occurs due to the hyperfine 
interaction neglected in our heavy-quark approximation. For simplicity we have considered only one D s state with 
mass 2.0 GeV. With the obtained potential matrix elements we solve numerically the eigenvalue equation for AE. 
It has essentially the same form as Eq. (|33[L but has now two terms due to the contributions from the two coupled 
channels. 

We find that the bare mass is shifted downwards by 51 MeV. Furthermore, the contributions from the various 
inelastic channels are given by 

AE(KD) = -35 MeV, AE(rjD s ) = -16 MeV. 

The magnitude of the total mass shift is clearly smaller here than found in the first considered case of the D*(2317) 
meson as can be explained by the presence of P-wave, rather than S- wave for the case of I?*(2317). It is clearly not 
sufficient to explain the experimental observed mass. We can also calculate the width, being the discontinuity of the 
right hand side of Eq. (33). We find T = 174 MeV, decaying predominantly into the K D-channel. The corresponding 
partial decay widths to the various channels are found to be 

T(KD) = 139 MeV, T(r]D s ) = 35 MeV. 

Clearly the predicted width is considerably larger than found in the SELEX experiment. We have also solved 
Eq. (33) for the case of the flavour symmetric value fx = f v = .093 GeV, which yields again a large value for the 
width. The above results suggest in this case that the resonance becomes very broad and does not support the SELEX 
observation, also the ratio of KD to r\D s channels quantitatively disagrees with the experiment. 
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In general, the size of the mass shift clearly depends on structure of the quark wave functions and hence it should 
be expected to be model dependent. Our study demonstrates, that the size of the shift due to channel coupling is 
in general large, but that it also can lead to very large widths in case that the resonance is above threshold of the 
coupled channels. As a result it can accommodate for the discrepancy between the predictions of dynamical quark 
models and the observed D*(1237) resonance, but there may exist situations where the resonance can become very 
broad due to inelastic channel coupling, as it is in the case of radial excited D* . 

V. DISCUSSION 

Let us compare our results to the existing in literature. The explanation of _D*(2317) as the cs p-wave level shifted 
down by the coupling to the decay channel D + K was considered in a series of papers [19], where authors have used 
a simple phenomenological model similar to our Eq. (26), (27) to describe the p-wave mesons. 

Another type of phenomenological model for the channel coupling, namely the model of Eichten at al. [26] was 
used in Ref. [27| to calculate the shift of the D s (0 + ) level and it was shown that the desired mass shift is obtained 
for a reasonable choice of parameters. 

Our results obtained with the parameter free chiral Lagrangian containing full x-dependence, qualitatively agree 
with those in Rcfs. [19] and [27] and exactly reproduce the experimentally found mass D* s (2318). 

We now turn to the state L>*(2632), found in Ref. ^H, but not yet confirmed by other groups (2^]. The theoretical 
prediction for the 2 3 S X state vary from 2774 MeV in Ref. [lj| to 2737 MeV made using the relativistic Salpeter 
equation in Ref. [2£j and 2716 in Ref. and in p rinciple are subject to the correction due to the global string 
breaking effect occurring for states of large size |3J| (m.s.r. radius of 2 3 Si state in Ref. 29] is around 1 fm). The 
expected correction is around - 20 MeV, which brings the theoretical mass of the D s (2 3 S'i) state to 2700-2720 MeV. 

Our calculation for the initial mass D*(l _ ), m = 2710 MeV, using Eq. (33), yields the shift AE = -76 MeV, 
with the total width T = 131 MeV, Ykd = HI MeV, Td sV = 20 MeV. Hence also in this case the total width and 
width ratio contradicts experimental data, implying that D* (2632) cannot be explained as the shifted l - level. 
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